Out-of-equilibrium heating of electron liquid: fermionic and bosonic temperatures 
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We investigate out-of-the equilibrium properties of the electron liquid in a two-dimensional dis- 
ordered superconductor subject to the electric bias and temperature gradient. We calculate kinetic 
coefficients and Nyquist noise, and find that they are characterized by distinct effective temper- 
atures: Te, characterizing single-particle excitations, Tcp, describing the Cooper pairs, and T^h, 
corresponding to electron-hole or dipole excitations. Varying the ratio between the electric j and 
thermal jth currents and boundary conditions one can heat different kind of excitations tuning their 
corresponding temperatures. We propose the experiment to determine these effective temperatures. 

PACS numbers: 74.45.-l-c, 73.23.-b, 74.78.Fk, 74.50.+r 
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A temperature, the term quantifying the common idea 
of "hot" and "cold," is one of the most fundamental con- 
cepts in physics. In statistical physics, the temperature is 
a parameter that controls the probability of the energy 
distribution of a given system over the possible states. 
The basic property of the temperature is that in an equi- 
librium it is the constant all over the system involved due 
to the interactions between the constituent subsystems. 
In out of the equilibrium case, different subsystems can 
acquire different temperatures. Electrons and phonons 
which equilibrate within themselves, but can possess dif- 
ferent temperatures if weakly coupled, is an exemplary 
case of such a non-equilibrium situation [^ || . 

In this Letter we show that strong correlations bring 
more richness to the non-equilibrium physics. We con- 
struct a kinetic theory of electronic transport in a disor- 
dered two-dimensional (2D) superconductor and demon- 
strate that distinct electronic subsystems respond differ- 
ently to external drive. We derive transport kinetic co- 
efficients and find that they are controlled by distinctive 
effective fermionic and bosonic temperatures each corre- 
sponding to a particular electronic subsystems: T^ for the 
single-particle excitations, T^^ for Cooper pairs, and T^h 
for electron- hole pairs (dipoles). These temperatures are 
defined through the respective energy distribution func- 
tions, and the relations between them are determined by 
the character of the external drive, or more specifically, 
by the ratio of its components, the applied voltage bias 
and/or temperature gradient (see Fig. 1). 

We consider 2D disordered systems such that the 
Thouless energy corresponding to diffusion across the 
film Ej^ = D/cP, where D is the diffusion coefficient, 
d is the film thickness, well exceeds all the fermionic and 
bosonic temperatures. Hereafter we will use units where 
ks = fi — c — 1 {c is the light velocity). The bias- 
dependent subsystem temperatures manifest themselves 
most profoundly in the leading quantum corrections 0] 
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FIG. 1. Temperatures of electron- hole pairs {T^h), Cooper 
pairs (Tcp) and quasi-particles (Tc) vs. the parameter x = 
R2/[Ri + R2), characterizing the ratio between tunneling re- 
sistances -Ri between the substrate and superconducting film 
and R2, between the film and the gate. The substrate tem- 
perature is Ti and the gate temperature is T2, T2/T\ = 0.1, 
the voltage drop between them is Vg- (a) Vg = 0, where only 
thermal current is present; (b) eVc/Ti = 3, with Voj > jth- 
Generally, in the film driven far from equilibrium and for 
R\ ~ R2, different interrelations between fermionic and 
bosonic temperatures can be realized: Te ~ Tc-p -C Toh for 
ith > VgJ, while Te ~ Tcp > Teh for jth < VgJ- 



to the sheet conductance in the singlet, 5G^^ (Aronov- 
Altshuler correction) and in the Copper channel, 5G"'^ 
(Maki-Thomson correction) close to the transition: 
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and in the leading contribution to the in-plane Nyquist 



Sth{i^ = 0) w 2auTch- 
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Here e is the electron charge, r is the elastic scatter- 
ing time, T(j, is the dephasing time, Tgl is the Ginzburg- 
Landau time indicating how close the system is to the su- 
perconducting transition, and a^ is the Drude conductiv- 
ity. We consider disordered system where T^hicp) ^ t^^ ^ 
e^/icT^d) < 1, and Tgl < r^. 

The fermionic and bosonic temperatures are defined 
through the respective distribution functions, that are 
assumed to be space and time independent. We start 
with the fermionic temperature To describing the low 
lying single-particle excitations. In an equilibrium, the 
electron (hole) gauge invariant distribution function in 



?(oq)/ 



/p(e) 



the disordered electron liquid is flfyiif) 
l/[exp(e/T) -I- 1], where e is the quasiparticle energy, 
and the temperature appears as its low-energy "cut-off" : 
T — e/2[l — 2f^f?i{e)]\^^Q. Generalizing onto an out of 
the equilibrium case, one arrives at: 



re = {-41ima,A(e)}- 

e->0 



(4) 



where ^(e) = [/c(e) -I- /h(e)]/2, and /o(h)(e) is a general 
nonequilibrium electron (hole) distribution function. 

The electron-hole (e-/i) pair temperature Teh emerges 
in the bosonic distribution function of the electron-hole 
excitations mediating, for example, the electronic trans- 
port in metals [H . The bosonic nature of the e-h envi- 
ronment is revealed by the fluctuation-dissipation the- 
orem (FDT) for the density-density correlation func- 
tion in a quasi-equilibrium regime: n^(q, oj) ~ 2i(l + 
2N^'^^) Imn^(w, q), where N^'^^ is the Bose distribution 
function of the e-h excitations, Ji^{x — y) is the retarded 
correlator of the electron density, and 11^ (x — y) is the 
Keldysh density correlator. The low-energy cut-off of the 
distribution function A^'^'i' defines Teh = w-^'i"'^ |w-!-o- In 
the nonequilibrium case we find 



Toh = lim ojN^, 



N.^, 



^ E j[^~h{e-)]Ue+)de, 

a—e.h 



(5) 
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where e± = eztuj/2. The integrand in (0) gives the prob- 
ability for the state with the energy e+ be occupied and 
the state with the energy e_ be empty. In an equilibrium 

Turning to the nonequilibrium generalization of the 
FDT in the Cooper channel, i.e. taking the Keldysh 
component of the fluctuation propagators in the Cooper 
channel, one arrives at the temperature for the fluctuat- 
ing Cooper pairs as 



-^cp — 21qZcp 



{/c(0)[l-/h(0)]+/h(0)[l-/e(0)]}. (7) 



Here zcp = (1 + (Te/fl)^)-! and ^-^ = 4£[/:,(0) - 
/R(e)]/7re^ where /n(e) = /e(-e) - /h(e). In an equi- 
librium Te = T, Zcp = 1, /f(0)[1 - /f(0)] = 1/4, and 
therefore Tcp = T. The combination of distribution func- 
tions appearing in the definition of Tcp can be viewed as 
the probability for generating a Cooper pair. 

We consider a setup shown in Fig. H, where the su- 
perconducting film is sandwiched between the substrate 
and the gate and is separated from them by the tun- 
nel barriers with the resistances i?i and i?2 , respectively. 
We can manipulate the density matrix of the system 
applying the gate bias, Vq, and changing the temper- 
atures Ti and Tj of the substrate and the gate. The 
electronic distribution functions defining the subsystem 
temperatures in the film are found from the respective 
kinetic equations. The thickness of the superconduct- 
ing film d is much smaller than the inelastic quasipar- 
ticle scattering length, thus the corresponding collision 
integral can be omitted. Then the kinetic equations 
in the stationary situation assume the form of the cur- 
rent conservation law, Vje(h)(e) = 0, where the spec- 



tral quasiparticle current je(h)(£) 



-DVf, 



(h)- 



Here 



we use the gauge where A = 0. The boundary condi- 
tions require the conservation of je(h) (e) and imply that 

Jc(h)(e) = [/ifh)(e) - fcll)i^)]/i'e-'^Rii where the indexes s 
and / denote substrate and film, respectively, and ly is the 
single-particle density of states at the Fermi level per spin 
projection H. Similar relation holds for je(h)(e) between 
the film and the gate. The gauge invariant distribution 
functions satisfy the kinetic equations which contain di- 
rectly the electromagnetic fields rather then their respec- 
tive potentials. In our case electric field enters the ki- 
netic equation through the boundary conditions. In order 
to obtain the gauge-invariant distribution functions, one 
uses the transformation /e(h)(e) = /e(h)(e i e(/)f), where 
4>f is the electric potential in the film, see, e.g., g, 0. 
Thus, 

/o(h)(e) = xfp{e±eVG{l -x),Ti) + 

il-x)MeTeVaX,T2), (8) 

where the temperature appears as the second argument 
in the Fermi-functions, and x = R2/RS with R^ = 
i?i -l- i?2- Expressing subsystem temperatures through 
the contact parameters and using definitions Eqs.(H)-(H), 
one finds: 



Teh ~ X^Ti + (1 - xfT2 + 

x{l — x)eVG coth 



eVr. 



yWTTrfy 



Te = 



X (1 — a;) 

^rich^ii^fi^^T2ch2^ 



(9) 
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The temperature Tcp is given by Eq. (R) , where Q, ^ be- 
comes 2Im{T2j:*'[l/2 - ieVdl - a;)/(27rTi)] + Ti(l - 



x)*'[l/2 + ieVGx/(27rT2)]}/(Tir27r2), and * is the 
digamma function. The Ginzburg-Landau time is de- 
fined as t-,i = 8zcpT,^k{e) - LM - U^,T,) + 
l/2]/e7r + 4zcpT2(l -2/H(0))/f7. At zero vohage, t-^ = 
8Te[a;ln(Ti/T2) — hi(Tc/T2)]/7r, Tc is the equihbriuni su- 
perconducting transition temperature. 

The net electric current has a standard form j = 
6^/e[Jc(£) — Jh(e)]/2. To derive the thermal current, we 
construct the gauge invariant vector out of j, ((/>, A) and 
the energy current j^; = v j^e[]e{e) + j/i(e)]/2: jth = 
J-E ~ <Ai- This agrees with the local equilibrium temper- 
ature gradient form of jth and satisfies the conservation 
law, Vjth = £h where £ is the electric field. In a thin 
film, the gradient is reduced to the difference at the in- 
terfaces, and in the setting of Fig.|| 



Jth 



2i?se2 



^(r,2-r|)-eV|(2x-i; 



(11) 



The same expression for the thermal current holds for 
the 2D film or wire connected through the edges and/or 
ends. Then the parameter x is nothing but the coordinate 
along the layer (wire) normalized by its width (length). 

The behavior of subsystem temperatures as functions 
of the external drive given by Eqs. ([7|)-(p^) is shown in 
Fig.|l|. Below, we concentrate on the case i?i — R2 
(i.e. X — 1/2), since then the system is the most strongly 
driven out of equilibrium. (Note that the equilibrium sit- 
uation is realized for x = and x = 1, with temperature 
being T2 and Ti, respectively.) If transport is dominated 
by the thermal current due to the temperature difference 
between the substrate and the gate, and t~^ <^T2 <^Ti, 
at zero bias Vo = 0, then at x = 1/2 the quasiparticles 
and fluctuating Cooper pairs are much colder than the 
electron-hole pairs: T^ w 2r2 ic Tcp < Teh « Ti/4. 
On the contrary, if the external drive is provided mostly 
by the voltage drop, the electron-hole pairs appear to 
become the coldest electronic subsystem for x = 1/2: 
To ~ Tcp > Toh, (see Fig. |]b). It reflects the fact 
that the electric current affects mostly charged parti- 
cles, rather than the neutral electron-hole pairs. The 
relation Tcp — Tc (see Fig.ll|a) holds at all x as long as 
Vg < min(ri, T2) [since zcp is of order of unity in Eq. (0) 
at small voltages]. The single particle excitations may 
become much hotter than the Cooper pairs, Tc ;» Tcp, 
at elevated Vq and away from x = 1/2, where Zcp signif- 
icantly deviates from unity, (see Fig.p^) g. 

We derive the kinetic coefficients using the Keldysh 
functional integral technique. The Keldysh parti- 
tion function in the coherent state basis is defined 
as: Z = J DtpDtljexp{iS['tp,^p]}, where 5[?/',?A] — 
Ic dt {U^c.^dt^c.) - H[^, V-]}. Here H ^ Ho + Hi^t is 
the Hamiltonian, C is the Keldysh contour and a =t)i 
is the spin variable. The one-particle Hamiltonian Hq = 
— /r i'a [(V — ieA)^ /2m + fJ. — Udis — ecj)] V'a, where A, 
(j) and Udis are vector, scalar and disorder potentials. 



and /J, is the Fermi energy; the tensor summation over 
the spin indices is implied. The interaction Hamilto- 
nian describes the electron-electron interaction in the sin- 
glet and Cooper channels, Hiat = (1/2) J^. ^, {pr^fl^' Pr' ^ 
(A/i^)(5r-r''/'t'/'iV'4.V't}j where F^^^ is the singlet channel 
interaction amplitude, p{r) = ipa{Y)'4'a{j^) is the local 
electron density. Averaging over Gaussian disorder and 
carrying out the standard decoupling in the four-fermion 
terms in the action via the Stratonovich-Hubbard fields 
{A,Q,$} ||l^, and integrating out the degrees of free- 
dom with the energies higher than l/r, we arrive at the 
Keldysh nonlinear cr-model action: 



S[Qic,^ic,^ic]^Sa + S^ + Sc, 



(12) 



Here 5a 



-Ki 



\)Tt[A]^YAic]/2X, S^ 



e'uTv[bKYh^]/2, Sq = iT:i^TT[D{drQKr/4. - EdtQK + 
ieb/cQic + lA/cQic] and Q^ = 1. We used the unitary 



limit F^_j., — !• 00 implying that the charge screening 
length in the electron liquid is much smaller than the 
mean free path. The action (^2|) holds while the fermionic 
and bosonic temperatures that follow from it are much 
smaller than l/r. The check mark above the field vari- 
able indicates that it is defined on the tensor product of 
the Keldysh and Nambu spaces spanned by the Pauli 
matrices ai and t^, i £ {0,x,y,z}, respectively. So, 
Y = (7x (^ tq, S — a-Q (g) Tz and S = $ -I- 0. Multi- 
plication in time space is implicitly assumed, and "Tr" 
includes an integration over real space. The subscript K, 
denotes the gauge transformed fields: $k: = '^ — dtJC, and 
Ak; = A+V^, where /C = [k'''(Jo+k'^ax]'i^To [A and $ are 
defined similarly]. Then A = [A'^'ctq -t- A'^a^] (E)f+- h.c, 
and Ajcir,t) = e*«2K:(r,t)^g-«eH/c(r,t)^ g^ ig defined in 
the same way. The quantum (q) and classical (cl) com- 
ponents are defined in the standard way as the half- 
sum and half-difference of the field values at the lower 
and upper brunches of the Keldysh time-contour. The 
field A becomes the superconducting order parameter 
on the mean-field level, while the saddle-point equa- 
tion for Q corresponds to the Usadel quasiclassical equa- 
tions where Q plays the role of the quasiclassical Greens 
function. The covariant spatial derivative is OtQk = 
'^vQk -ie[SAK;,QK;]- 

We use linear response formalism to find the interac- 
tion corrections to the kinetic coefficients of the electron 
liquid integrating out the fluctuations around the metal- 
lic saddle point of S", A = U{(Jz ® fz)^~^, where 



Ut.t'{v)=U-l,{v) 







-^^t-t'+QTo 



(13) 



Here Ti_t'(r) = d\a,g{F^ ^, {r) , Fl\, {r)) . Performing the 

Wigner transformation we map F^\, to the quasiparti- 
cle distribution functions: T'o(h)(e,T) = 1 — 2/e(/j)(e, r). 
We split the gauge field, k, and electromagnetic fields into 



the slow, (k), and fast, Sk, components. The fast com- 
ponents describe the fluctuations in the electron system 
and slow components are related to the gauge transfor- 
mations of the external field potentials M . To optimize 
the fluctuations of the quasiparticle phases induced by 
the Coulomb forces ||, |l^ we solve the equation that 
couples Sk with the fluctuating electromagnetic fields: 
Sbl - DV Sg'i = 0, 

DV5g'\r, Lo) + Sb'^ir, lo) - -2(1 + 27V„)DV,5g«(r, uj) , 

(14) 



where 5b°^''^^ = ,5$<i('=i) - dt5k'i'^''^\ (5g<i('=i) = V(5/fci('='). 
This way we fix the gauge and define N^ as given by 
Eq. (g). In the local equilibrium, /'^'^'''(e, r) = /j^(e T 
e{(j) — 9r (fc'^'))), and then choosing (fc'^') = (j)t one obtains 
the gauge invariant distribution functions. 

Now, the electric conductivity is obtained by differ- 
entiating the partition function over the quantum and 
classical components of the vector potential, and then in- 
tegrating out the gauge. A- and Q-fluctuations (diffusion 
and Cooperon degrees of freedom). The major fluctua- 
tion contribution to the in-plane conductivity close to the 
transition is the Maki-Thompson correction: 



Sa 



MT 



27r 



[L-']kX 



q,cj,€ 

2 



LR,Uq)rdJ,ie) + |Lfl.,_^(g)p9Jh(£) 



(15) 



where [L ^]k ~ iT:Tcp/{2TcZcp), Lj^^^ = 7r[-(TGLZCp) ^ + 
{^Dq'+iLo-2ie(f>f-iATeil-2U{0))){l+iT,n-^)]/{8T,). 
In two dimensions the integral (Rq) diverges loga- 
rithmically at small q, and should be cut off at 
l/r^ leading to the result (|^), for Tgl <C t^ = 
aodn^el + T-:') / {8e^Tc,TJ{eo)). Here /(eq) - 
(7r|eo| - 2eoarctan[r-^ie^^] + t-^ log[4T-LV (eg + t-;^)]), 
where eo = 4Te(l - 2/^(0)) - T~^Te/il.. Corrections from 
the fluctuating Cooper pairs are suppressed by the mag- 
netic field perpendicular to the film, and the Aronov- 
Altshuler correction becomes the dominating one: 



Sa w Im 



q.o; 



d^ Hi + 2JV^)] 



(16) 



This integral also diverges logarithmically in 2D, and the 
integration over uj should be cut off at the upper limit by 
1/t pT[ , whereas the infrared boundary is determined by 
Teh g, giving rise to Eq. (0)|l|. 

Next, we calculate the thermal Nyquist noise. In the 
noninteracting case, the current density j = cTniC + ^]- 



Here ^ is the random electrical field induced by the fluc- 
tuations of the vector potential quantum components, 
the statistical properties of which are defined by the 
Fourier transform (C^)tj = '^^(l + 2N^)/a-o. Then, we 
obtain the Nyquist noise S'th(i^) — crDw(l + 2Ni^). This 
result is derived in the limit V -^ u. If w <C Teh then 
[1 + 2N^] -)■ 2Teh/w, and we get Eq. ^. 

Fermionic and bosonic temperatures can be straight- 
forwardly detected, by direct measurements of the mag- 
netoresistance and temperature dependence of the resis- 
tance at fixed magnetic fields using the setup shown in 
Fig. 1. Varying Ti and T2 and the gate voltage, one can 
directly observe the effect of Tcp in the Maki-Thompson 
correction. At the same time changing temperatures Ti 
and T2 at the fixed high magnetic field, where the Maki- 
Thompson contribution is suppressed, one infers the in- 
formation about Teh. Also T^h can be determined from 
the Nyquist noise measurements. 
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